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AIRFOIL PROFILES FOR MINIMUM PRESSURE DRAG AT SUPERSONIC VELOCITIES— GENERAL 
ANALYSIS WITH APPLICATION TO LINEARIZED SUPERSONIC FLOW ^ 

By Dean R. Chapman 


SUMMARY 

A theoretical investigation is made of the airfoil profile for 
minimum pressure drag at zero lift in supersonic flow. In the 
first part of the report a general method is developed for calcu- 
lating the profile having the least pressure drag for a given 
auxiliary condition, such as a given structural reguirement or a 
given thickness ratio. The various structural requirements con- 
sidered include bending strength, bending stiffness, torsional 
strength, and torsional stiffness. No assumption is made 
regarding the trailing-edge thickness; the optimum value is 
determined in the calculations as a function of the base pressure. 

To illustrate the general method, the optimum airfoil, defined 
as the airfoil having minimum pressure drag for a given auxil- 
iary condition, is calculated in a second part of the report using 
the equations of linearized supersonic flow. It is found that the 
optimum aiifoil in most cases has a blunt trailing edge. It also 
is found that the optimum thickness distribution dejmnds only on 
one dimensionless parameter, termed the “base pressure param- 
eter” . This parameter involves the Mach number, airfoil 
thickness ratio, and base pressure coefficient. The effect of 
variations in each of these latter three quantities on the shape of 
the optimum profile is discussed _ and a simple critenon formu- 
lated for determining the condition under which the optimum 
trailing-edge thickness is greater than zero. The calculated 
pressure drag of the optimum profile is compared to that of a 
biconvex sharp-trailing-edge profile satisfying the same struc- 
tural requirement. The reduction in pressure drag depends on 
the base pressure parameter and varies from a few percent to as 
mibch as 75 percent. 

INTRODUCTION 

In supersonic flow the finite thickness of an airfoil invari- 
ably introduces a certain amount of pressure drag which can 
be minimized by a rational choice of airfoil shape. The 
profile for minimum pressure drag depends, among other 
things, on the particular auxiliary condition that is imposed 
on the airfoil geometry. For example, if it is required that 
the optimum profile (defined herein as the profile of least 
pressure drag for a given auxiliary condition) satisfy the 
auxiliary condition of a given thickness ratio, then according 
to a well-known result of Ackeret’s linearized airfoil theory, 
the so-called double-wedge profile represents the optimum 


sharp-trailing-edge airfoil. This particular auxiliary condi- 
tion, however, does not represent practical cases where an 
airfoil must satisfy a certain structural requirement, such as 
a given area moment of inertia, or a given section modulus. 
Drougge (reference 1) has made a more elaborate theoretical 
analysis to determine the optimum profile for the auxiliary 
condition of a given bending stiffness of the airfoil, and also 
for the condition of a given torsional stiffness. Drougge used 
linearized airfoil theory and considered only sharp-trailing- 
edge airfoils. His results are somewhat limited in two 
respects: They do not cover cases outside the scope of linear- 
ized airfoil theoiy, and, though they include the au.xiliaz’y 
conditions of given bending and torsional stiffness, thc 3 '^ do 
not include the auxiliary condition of a given bending strength 
(given section modulus). A far more important limitation 
of this analysis, though, is the tacit assumption that the 
optimum airfoil will have a sharp trailing edge. 

There is a small amount of experimental evidence in the 
measurements of Ferri (reference 2) on airfoils witli shaip 
trailing edges which suggests that the optimum airfoil might, 
in fact, have a moderately thick trailing edge. The meas- 
ured profile drag of one airfoil tested by Ferri (G. U. 3 airfoil 
at a Mach number of 1 .85) was considerably lower than 
inviscid theory would indicate. Allowance for skin friction 
would cause this discrepancy to become even greater. 
Schlieren photographs and pressure-distribution measure- 
ments showed that viscous effects effectively thickened the 
airfoil shape near the trailing edge. From these results it 
can be inferred that at moderate supersonic velocities it is 
possible for an airfoil with a thickened trailing edge — that is, 
a blunt-trailing-edge au'foil — to have lower drag than a 
corresponding sharp-trailing-edge airfoil. Employing a 
different approach, this inference has been obtained from 
quantitative considerations in reference 3, where a reasonable 
estimate of the base pressure was made and the drag calcu- 
lated as a function of trailing-edge thiclaiess. Such calcula- 
tions, though very approximate, have indicated that ui cer- 
tain cases a moderate increase in trailing-edge thiclaiess will 
decrease the over-all pressure drag. 

Apart from the reasons just cited for expecting that the 
optimum supersonic airfoil might have a thick trailing edge, 
there are other independent considerations which suggest 


* Supersedes NACA TN 22G4, “Airfoil Profiles for Mitiimum Pressure Drag at Supersonic Velocities— General Analysis With Application to Linearized Supersonic Flow,” by Dean R. 
Chapman, 1951. Various examples of optimum profiles given in TN 2284 have been supplemented and revised for the present report in accordance with experimental measurements of base 
pressure published subsequent to TN 2264. 
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the same result. By analyzing conditions at infinite Mach 
number, Saenger pointed out in 1933 that even with a vacuiun 
at the base the optimum airfoil for a given thickness ratio 
would, in this extreme case, have a traUing-edge thickness 
equal to the maximum airfoil thickness. (See reference 4.) 
In reference 5, Ivey obtained a similar result by calculating 
the pressure drag at a Mach number of 8 for a family of 
airfoils having various positions of maximum thickness. 
More recently, Smelt (reference 6) has developed an approx- 
imate condition determining when an au-foil with maximum 
thickness at the trailing edge has lower drag in hypersonic 
flow than an airfoil with a sharp trailing edge. Saenger, 
Ivey, and Smelt, however, did not consider airfoils having a 
trailing-edge thickness less than the maximum airfoil thick- 
ness, and hence their results do not determine the optimum 
profile for hypersonic velocities. Nevertheless, it is evident 
that at high supersonic Mach numbers the optimum profile 
has a relatively thick trailing edge. On this basis it is not 
unreasonable to expect that at lower supersonic Mach num- 
bers the optimum profile would have some thickness at the 
trailing edge. 

The physical reason why it is possible for a blunt-trailing- 
edge airfoil in supersonic flow to have a lower pressm’e drag 
than a corresponding sharp-trailing-edge airfoil is quite 
simple, as can be illustrated by the two profiles shown in 
figm’e 1. These profiles have the same area, which corre- 
sponds to the same torsional stiflfness of a thin-sldn structure 

/ Sharp- traiUng-ec/ge ait 'foil 


r - - Blunf-frailing-edge airfoil 



Figure 1, — Sketch comparing a typical sharp-trailing-edge airfoil and a blimt-trailing-edge 
airfoil of equal area (equal torsional stiffness for a thin-skin structure). 


The blunt- trailing-edge airfoil has a slightly smaller thick- 
ness ratio and a position of maximum thickness which is 
farther rearward, hence the leading-edge angle is smaller 
and the pressure drag of the surface forward of the trailing 
edge is less than that of the sharp-trailing-edge airfoil. A 
certain amount of base drag, however, obviously is added by 
employing a thick trailing edge. If the added base drag is 
less than the reduction in pressure foredrag, then the net 
result is a smaller total pressure drag for the blunt- trading- 
edge airfoil. This invariably is the case at extremely high 
supersonic Mach numbers where the base drag is negligible 
compared to the pressure foredrag. At low supersonic 
Mach numbers, though, the base drag can be many times 
the pressure foredrag, and the optimum trading-edge thick- 
ness must be expected a priori to depend to a great extent 
on the base pressure. 

The present theoretical analysis was initiated in view of 
the foregoing considerations. The primary purpose of the 
investigation is to develop a method of determining the 
supersonic airfoil profile for minimum pressure drag at zero 
lift, without making an arbitrary assumption about the trail- 


ing-edge thickness. The profile so determined, which is 
termed an optimtun profile, is considered to depend on the 
base pressure, Mach number, and the particular auxiliary 
condition imposed on the airfoil. A secondary purpose of 
the investigation is to develop a method of sufficient gener- 
ality to enable second-order and shock expansion theories 
to be used in calculating optimum profiles. Such generality 
is desirable in order to obtain results that are valid at hyper- 
sonic Mach numbers. 


NOTATION 


B base pressure parameter for linearized supersonic flow 

B limiting value of the base pressure parameter below 
which the optimum airfoil has a blunt trailing edge 
c airfoil chord 

Ca section pressure drag coefficient 

/ symbol for the function Py' 

F(k, <p) incomplete elliptic integral of the first kind of modu- 
lus k and amplitude <p 
h trailing-edge thickness 

H dimensionless trailing-edge thickness {hjt) 

I given value of the auxiliary integral [( r wf 

kn constant defined by equation (20) 

I length of chord over which airfoil thickness is constant 
L dimensionless length of chord over wliich airfoil thick- 
ness is constant (Z/s) 

M Mach number 

n arbitrary parameter appearing in the definition of the 
auxiliary integral I 

(For the examples considered n is taken as 1, 2, 3 or 
”•) 

p static pressure on airfoil surface 

P pressure coefficient j^(p— p„) j Pc„ 

base pressure coefficient 

Be. Reynolds number based on airfoil chord 
s distance from leading edge to first downstream posi- 
tion of maximum thickness 
t maximum thickness of airfoil 

V velocity 

X distance from leading edge 

X dimensionless distance from leading edge (x/s) 
y ordinate of upper half of airfoil 

Y dimensionless ordinate [y/{t/2)] 

^MJ~1 

X Langrangian multiplier (arbitrary constant) 

(T arbitrary parameter appearing in the definition of the 
auxiliary integral I 

(For the examples considered a is taken as 0 or 1.) 
p mass density 

SUBSCRIPTS 

0 airfoil surface at leading edge 

1 airfoil surface at trailing edge 

CO free stream 

h base 
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ca circular-arc biconvex airfoil with sharp trailing edge 
d/m double-wedge airfoil with sharp trailing edge 

V vacuum at base 

SUPERSCRIPT 

' differentiation with respect to x 

THEORETICAL ANALYSIS 
ASSUMPTIONS AND STATEMENT OF PROBLEM 

In the analysis which follows several simplifying assump- 
tions are made. Two-dimensional airfoils in a purely super- 
sonic flow at zero lift only are considered. It is assumed that 
the pressure at any point on the airfoil surface forward of the 
trailing edge can be calculated from the flow of an inviscid, 
nonconducting gas. It is further assumed that the leading 
edge is sharp. No analogous assumption is made regarding 
the trailing-edge thickness, but it is assumed that the base 
pressure coefficient P\, is known. This enables the optimum 
trailing-edge thickness to be calculated as a function of Ps; 
hence, experimental data on base pressure in two-dimensional 
flow are required in order to apply the theoretical results of 
the analysis to a given case. 

From the fact that the surface pressures on the top and 
bottom of an airfoil can be calculated independently in a 
supersonic flow, it follows that at zero lift the optimum 
profile will be symmetrical about the chord plane. Conse- 
quently, reference is made throughout to the thickness distri- 
bution of only the upper surface of an optimum profile. 

In comparing the pressure drag of various profiles, the 
chord length is held constant, and the thickness distribution 
along the chord is varied in a manner which is arbitrary 
except for the requirement of satisfying the particular 
auxiliary condition being considered. The various auxiliary 
conditions investigated are: a given torsional stiffness of the 
airfoil section, a given torsional strength, a given bending 
stiffness, a given bending strength, and a given thickness 
ratio. For each of the structural conditions the case of a 
thin-skin structure and a solid-section structure is considered 
since the optimum airfoil profile may be expected to depend 
somewhat on the type of structure. Attention is focused on 
the fact that the basic idea employed in the analysis involves 
the minimizing of pressure drag for a given structural require- 
ment; the results obtained with this method of approach are 
the same as would be obtained if the structural characteristic 
were maximized for a given value of the drag.“ 

MATHEMATICAL FORMULATION OP PROBLEM 

The pressure drag of an airfoil with a thick trailing edge 
is the sum of the base drag and the pressure drag of the sur- 
face forward of the trailing edge. Letting P be the surface 
pressure coefiicient, y{x) the function defining the surface, 
and Pj the base pressure coeffieient, then may be expressed 
as 

fPy'dr-Pj- (1) 

C Jo c 

The problem is to determine the particular function y{x) 
and the corresponding value of the trailing-edge thickness 


h which minimizes this expression for a given auxiliary 
condition. 

Before expressing the various auxiliary conditions in ana- 
lytical form, it should be noted that the surface pressure co- 
eflScient P is regarded as a known function of the variable y' 
and the two parameters y\ (sm-face slope at leading edge) 
and The actual functional form of P{y',y'o,M„) will 

depend on whether linearized, second-order, hypersonic, or 
shock-expansion theory is employed in calculating surface 
pressmes. For example, if li near theo ry were employed, the 
explicit expression P=2y ' — 1 would be used; but, 
if shock-expansion theory were employed, a more complex 
implicit expression involving y'o as well as y' and Ma> would 
have to be used. In order to allow various theories to be 
employed, the particular functional form of P{y' ,y\,MP) 
will at present be unspecified. The equations which result 
can be applied to any of the various theories by substituting 
the appropriate function for P. 

Turning now to the consideration of auxiliary conditions, 
it is clear that some integral expression will be involved, 
since the function y{x) is not known beforehand. If, for ex- 
ample, the airfoil is a solid-section structure and the moment 
of inertia is prescribed, then the particular auxiliary condi- 
tion which y{x) must satisfy in addition to minimizing Ca is 

that the integral y^dx be constant. A different auxiliary 

condition would, of course, be represented by a different 
integral. In the present investigation a somewhat general- 
ized auxiliary condition is used which is represented by the 
single integral 

I=— 77 ^ da: = constant (2) 

cjo {tj2y 

where n and v are constants. Thus the example just cited 
is a special case of the above integral with n=3 and o-=0. 
To illustrate further, the auxiliary condition of a given section 
modulus of a solid-section airfoil is represented by the special 
case n,= 3 and o-=l. The corresponding solution for y(x) 
in this latter case would provide the profile of least pressure 
drag for a given bending strength. 

Some of the different structural criteria to which the gen- 
eral integral (2) eorresponds are summarized in the following 
table:® 


n 

<r 

Structural criteria 

1 

0 

Given torsional stiffness, or torsional strength, or volume of thin- 



skin structure 

2 

0 

Given bending stiffness of thin-skin structure 

3 

0 

1 Given bending stiffness or given torsional stiffness of solid-section 



structure * 

2 

1 

Given bending strength of thin-skin structure 

3 

1 

Given bending strength of solid-section structure 


• As a first approximation the torsional stiffness of a thin solid*section profile is taken to 
be proportional to the moment of inertia about the chord plane. 


Thus, by solving the problem with the general integral 
(2) left in terms of n and a a wide variety of auxiliary con- 


2 This statement, which appears evident from physical considerations, is equivalent to Mayer’s reciprocity theorem for isoperimetric problems in the calculus of variations. 

3 For thin-skin structures the thickness of skin is taken to be constant over the chord length. The two cases ti=o-= 0 and 7i=a-=l are not included in this table as they apparently represent 
no sensible practical problem. 
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ditions can be obtained simply by substituting appropriate 
integers for n and <r. From an engineering viewpoint the 
general form of equation (2) enables approximate solutions 
to be obtained for wings of intermediate structural solidity 
by interpolating between the solution for essentially zero 
solidity (thin-skin structures) and the solution for complete 
solidity (solid-section structure). 

Summarizing, the problem formulated can be stated math- 
ematically as that of finding the airfoil-ordinate function 
y{x), and the trailing-edge thickness h—y(c), which mini- 
mizes the drag expression (1) for a given constant value of 
the structural integral (2). The boundarj^ conditions im- 
posed are that y{Q) = 0 and that Pt, is given. If t did not 
appear in equation (2), this mathematical problem would be 
a relatively simple isoperimetric problem in the calculus 
of variations. The occurrence of t, the maximum value of 
y(x), complicates matters because it is not known before- 
hand and, in fact, is one of the quantities to be determined 
from the given values of M„, Pj, and I. Actually, all 
equations necessary for solving the problem formulated could 
be obtained directl}’ from advanced treatises on the calculus 
of variations since the problem is a special case of the so- 
called “problem of Bolza with variable end points”. (Such 
a procedure would lead quickly to equations (11), (12), and 
(13).) However, the necessary equations can also be ob- 
tained bj^ the simple methods employed here. 


Application of this basic principle, as will be seen, leads to 
a sufficient number of equations to determine the complete 
geometry of the optimum profile. 

A sketch of the type of optimum profile to be analyzed 
and the corresponding varied profile is shown in figure 2. 
Various quantities which appear often in the subsequent 
analysis are illustrated in this figure. It is to be noted that 
allowance is made for the possibility that the optimum 
profile may have a straight midsection of length I, the opti- 
mmn value of which must be determined from the analysis. 
The varied profile is selected such that it does not change 
the ordinate or the surface slope at the leading edge. Intro- 
ducing the definition Py' for the purpose of brevity, and 
equating the drag of the optimum profile to the drag of 
the infinitesimally varied profile, yields 




h-\-Sh 

b — o — 


The small change in slope 8y' is equal to ^ (Sy), so equa- 
tion (3) can be written 

Intergrating bj^ parts. 


' 5 / d(^ , a 

dx “ ^ " 2 


(4) 


METHOD OF SOLUTION 

Given structural criteria. — Since the pressme drag of the 
optimum airfoil, by definition, is the least possible of all 
airfoils having a given value of the structural integral (2), 
it follows that the pressure drag of any “varied” airfoil, 
having ordinates and slopes everywhere close to those of 
the optimum airfoil, must be in the neighborhood of a 
minimum. Hence, by considering only infinitesimal changes 
8y in the ordinate of the optimum profile, the corresponding 
increment in drag of such a varied profile can be equated 
to zero. Since y{x) is to provide the true minimum, the 
resulting equation must hold for an arbitrarv ordinate change 
5y varying with x, or for an arbitrary change in airfoil thick- 
ness 5t, or for an arbitrary change in trailing-edge thickness 
8h, or for any combination of variations thereof, provided 
only that the integral (2) is constant for all such variations. 


0 = 



X=C 

x—0 



dx 



dx — 




dx—Pu 


2 


or, finally 


5Crf=0 = 






A 

dx VbyV 




(5) 


The chordwise distribution of the variation 8y is not en- 
tirely arbitrary; it must be such that the auxiliary condition 
is satisfied, namely, the value of I for the optimum profile 
must be equal to that of the varied profile. 


i r y" 

C Jo(G2)' 


dr = /=- 
c 


p {y+8yY 
Jo(t/2-|-5</2)- 


dx 


(6) 


\ — Optimum profile 



Retaining only first order variations, this e.xprossion simplifies 
to 

J y"dx= j' y"dx + J" ny”~^ 8ydx — <7^^ y"- dx (7) 
or 

0 = (^0 j ( 8 ) 

This equation must be satisfied, of course if both terms on the 
right side are multiplied by an arbitrary constant X. More- 
over, equation (5) must be satisfied simultaneously. The 
arbitrarj' character of X enables the two equations (5) and (8) 
to be combined into a single equation which must hold for 
arbitrarj' variations in 5y, 8h, and 8t. 



AIKFOIL PROFILES FOR MINIMUM PRESSURE DRAG AT SUPERSONIC VELOCITIES 


5 




ah 


+ 


\acl 



( 9 ) 


From this it is seen that a first integral of the basic differential 
equation (11) is 

hP 

y'^ 7 + Xj/“= constant (1 la) 


Thus, since X is arbitrary, this single equation implies that 
both equations (5) and (8) are satisfied. 

If at this point the variation at is considered to be arbitrary, 
equation (9) as written would incorrectly suggest that 
\acl must be zero. Hence it is to be expected that an 
additional term containing 8t exists in the integral expression 
of equation (9). Such a term arises from the contribution 
of the straight midsection to this integral, since over this 

d /d/\ 

region, 5y= 6^/2. Also, y=i/2 and J=0 over this region. 

Hence 




The variations 8y, 8h, and 8t can now be conducted entirely 
independent of each other. Each of the bracketed terms in 
equation (10) must be zero, if the individual variations are 
not zero. Reemploying the definition /=Py', the following 
equations are obtained: 


for 8ij9^0 
for 8h9^0 
for \8t7^0 



(11) 


(12) 

1 I 

c (J/2)"-" ° 

(13) 


The differential equation (11), of course, results from 
equating to zero each of the two integrals in equation (10). 
This differential equation, therefore, need be satisfied only in 
the two chordwise regions covered by the limits of these 
integrals, namely, in the region from x=Q to x=s, and in the 
region from x=s-\-l to x=c. (See fig. 2.) If the optimum 
airfoil has a finite length of straight midsection (e.g., AB in 
fig. 2), the differential equation (11) need not be satisfied in 
this intermediate region. 

Fortunately, one integration of equation (11) can im- 
mediately be made, thereby lowering the order of the basic 
differential equation to be solved. Multiplying equation 
(11) by y' gives 




— 2v'v" ^ ^ 




At the point , or points, where y'=0 the ordinate is equal to 
t/2. Evaluating the constant of equation ((11a) from this 
consideration yields 

This equation, together with equation (12), equation (13), 
the given value of /, and the boundary condition i/(0)=0, 
determines the complete geometry of the optimum profile. 

Given thickness ratio. — Attention is called to the fact that 
special precautions must be taken in applying the foregoing 
analysis to the auxiliary condition of a given thickness ratio. 
For this particular case 8t is zero, thereby causing the last 
term in equation (10) to vanish automatically without 
requiring equation (13) to be satisfied; equation (13), there- 
fore, does not necessarily apply when the thickness ratio 
is prescribed. Moreover, equation (11b) also does not 
necessarily apply since it was assumed in the process of 
obtaining this latter equation that the optimum airfoil had 
at least one point where y' = 0. Such is not the case for 
the auxiliary condition of a given thickness ratio, and hence 
more detailed consideration is required. 

The appropriate differential equation to be used when 
tje is given may be obtained from equation (lla)^ by setting 
X=0. There results 


^^=constant (14) 

which is satisfied by any straight surface y' = constant, 
regardless of whether linearized, second-order, or shock- 
expansion theory is used for P. The appropriate condition 
which must be satisfied at the trailing edge is, from equation 
( 10 ), 


[(P+!/'^)-P.]**aO (16) 

Here the inequality is included since 8h for the case of a 
given thickness ratio is not always entirely arbitrary. Thus, 
when h—tl2 (wedge airfoil) 8h is restricted to always be neg- 
ative, and a minimum can exist if 

(f), 

this would make 8ca always positive instead of just 
making stationary. Consequently, under certain con- 
ditions two solutions are possible. First, the upper half 
of the optimum profile may consist of two straight segments 
with (as illustrated in fig. 3), provided the equal sign 
in (15) apphes. Second, the optimum profile may be a 
wedge profile with h=t, provided the inequality sign in (15) 
apphes. If both types of solution are physically possible 


< If equation (11) is used there results P-f-|f'dP/&y'=constant, which also is satisfied by any straight surface. The constant in this latter equation, however, does not have the same value 
for both straight segments comprising the profile; whereas, the constant in equation (14) is the same for both segments. (See appendix.) 
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Figure 3.— Sketch of upper half of optiraum profile for a given thickness ratio. 


same thiclcness as the optimum airfoil, it is evident that for 
an actual airfoil //(t/ 2 )”“‘' will be of the order of one-half or 
two-thirds. As an example, this means that when bending 
strength is critical in a thin-skin structure (n=2), the opti- 
mum length of straight midsection will be of the order of 
one-third the chord length. 

CALCULATION OF OPTIMUM PROFILES USING LINEAR 
AIRFOIL THEORY 


in a given case, the true solution, of comse, would be the 
one with lower drag. 

If it is possible to obtain a general solution to equations 
(11b), (12), and (13) (such is the case for linearized super- 
sonic flow), then the first of the two solutions mentioned 
above may be obtained without solving equations (14) and 
(15), but by passing the general solution to the limit as 
71— > 00 ,5 In order to verify that this limiting condition 
represents the auxiliary condition of a given thiclmess ratio, 
three conditions must be satisfied: First, as w— >oo the 
auxiliary integral must correspond to the case of a given 
thickness ratio; second, the differential equation (lib) must 
reduce to equation (14) in the limit as ; and third, 

the infinite value of n must be compatible with equation (13). 
Consider for the time being that the chord is of unit length. 
Since for any reasonable airfoil y<Cl, it follows that 


r 


y” dx -^0 as n-^ <» 


and hence 7-^0. 


A solution for 7—0 


would represent the optimum airfoil determined without regard 
for an auxiliary integral. Such is the condition that would 
be used in determining the optimum airfoil for a given 
thickness ratio; hence the first of the above-mentioned 
conditions is satisfied. Inasmuch as y/(</ 2 )<(l, it is evident 
that [t//(</ 2)]"— ^0 as thus reducing the differential 

equation (lib) to the form 


y 


n 


'6P 



—constant 


which is the correct differential equation. As regards the 
third condition to be satisfied, it is seen that the values 
7=0, 71 = 00 , (= 0 , and o-= finite are compatible with equation 
(13). Consequently, the limiting case -a— > in the general 

solution to equations (11b), (12), and (13) represents one of 
the possible solutions for the case of a given thickness ratio. 
T his fact will be used later in the report. 

QUALITATIVE RESULTS OBTAINABLE WITHOUT SPECIALIZING TO A GIVEN 
TYPE or SUPERSONIC FLOW 


SPECIALIZATION AND SOLUTION OF GENERAL EQUATIONS FOR LINERAR- 
IZED SUPERSONIC FLO W 

Given structural criteria. — According to the theory of 
linearized supersonic flow, the local surface pressure coeffi- 
cient on an airfoil is given by 

P=\y' 


For this approximation the basic differential equation (11b) 
becomes 


7/^| = X[(</2)''-7/«] 
or, after solving for dxjdy, 


V¥ 


x = 


r dy 


n 


(17) 


(18) 


This can be put into a more convenient form by introducing the 
dimensionless variables X,Y,L, and H defined in the lis t of S 5 ’’m- 
bols, and eliminating .^X /3 by evaluating equation (18) ata:=s. 
Between the leading edge and the first downstream position 
of maximmn thickness dy/dx is positive while the dimension- 
less variables X and Y both vary from 0 to 1 . Along the 
length of straight midsection Y—1, while X varies from 1 to 
l-f-A. Along the downstream portion of curved surface 
dyidx is negative and Y varies from 1 to H, while X varies 
from l-\~L to c/s. Consequently, equation (18) giving X as a 
function of Y becomes (with the convention that the sign of 
all radicals is positive) 




1 — ^ surface facing upstream 

- 7 ^==^ on surface facing downstream 


(19) 


Although few quantitive results can be obtained from the 
basic systems of equations (11), (12), and (13), without 
specifying a particular form for the surface pressure coefficient 
P, there is one general qualitative result that can be obtained 
from equation (13) without any further calculation. The 
optimum length of straight midsection I always is zero for the 
auxiliary condition of a given torsional stiffness or a given 
bending stifl'ness (o-= 0 ), but never is zero for the auxiliary 
condition of a given bending strength (o-=l). Since (</ 2 )’'“" 
would be the value of 7 for a rectangular-bar airfoil of the 


where the constant depends only on n and is given by the 
definite integral 

/2 for n= 1 

k„= f =< tt/ 2 for 77=2 (20) 

JoVl-F'‘ ) 

U.4023 . . . for ?i = 3 

It may be noted here that integrals of the type occurring 
in equation (19) also occur at numerous places in the sub- 
sequent analysis. Such an integral, being a function of the 


5 It should be noted that the value n=0, when substituted in the auxiliary integral, gives /==(</2)-*=constant; but this value cannot be used to obtain the solution for the case of a given 
thickness ratio because n-0 is incompatible with equation (13). 
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parameter n and the lower limit Y, can be evaluated either 
graphically or analytically. An analytical evaluation for 
the first three integer values of n yields 


and 


L= 


2(7 

k„(n—a) (»+2) 





where 



2-sjl—y for 71=1 
7r/2— arc sin Y for n=2 
</>) for 71=3 


ib=sin y|=0-9659 . 


cos ^0 = 


V3-1 + F 
V3 + 1-F 


(21) 


The function F {k, <p) is the incomplete elliptic integral of the 
first kind of modulus k and amplitude (A table of this 
function is given in reference 7, page 122.) For convenience 
the v^arious formulas developed later are left in terms of the 
above integral; specialization to the individual functions 
indicated in equation (21) for a given n could be made in 
any subsequent formula if desired. 

For linearized supersonic flow, equation (12) becomes 

p.=p.+(’/ ^1), 

/^\ 

^ \dx/i 
_2(i/c) {dY\ 

Ksic) \dX), 


or, on using the relation —k„ dX=dY!^jl —■ F“ which applies 
to the surface facing downstream (equation (9)), there is 
obtained 


P,<3 2kJl-H’‘ 
tic (s/c) 


( 22 ) 


As defined earlier, H = hlt is the optimum trailing-edge 
thickness expressed as a fraction of the maximum thickness. 
Equation (13) can be written, in terms of X and Y, as 


-=-( r Y’^dX+l+ C" Y'^dx) (23a) 

s V Jo Jl+i / 


or, after specializing to linearized flow. 



/ Y’^dY , j., f Y'^dY \ 


(23) 


Equations (22) and (23) can be put into more usable forms 
by noting from equation (19) that 



dY 

vr^ 


(24) 


After some algebraic manipulation involving integration by 
parts and introduction of the definition B = — Pi,;3/(t/c) 
there results from combining equations (22), (23), and (24), 


27i,(7i,-l-2 — (7)Vl — 

2<r27Vl-2?« 

ri dY 

{n—a) {n-\-2,) 

” ' 7l(7l + 2-cr) ' 

Ju Vl-F“J 


(25) 


These latter two equations are the final equations deter- 
mining the optimum dimensionless trailing-edge thickness 
and the optimum dimensionless length of straight midsection 
L. The corresponding equations involving the given value 
of I can be developed from equations (13), (22), (25), and 
(26) as follows: 


I In s / T n\ 2fc„Vl— / r 

?) 

dY 

Vl — F" 

(i/2)”"'*^ € <T c\ 

2 

-rr B 

•• 

71-\~2 — O' 

, 2.^71-27" 

p dY 

" ' n{n^2-<r) ' , 

Jy 7l_F" 




This last equation determines i? as a function of 7/(i/2)’‘“'', 
or vice versa. It is to be noted from equations (24), (26), 
and (27) that the geometry of an optimum profile for given 
values of n and a is determined solely by H, which, in tm’n, 
depends only on the base pressure parameter B (equation 
(25)). 

Given thickness ratio — Since the use of linerarized theory 
provides a general solution in closed form of the basic equa- 
tions (11b), (12), and (13), the optimum profile for a given 
thickness ratio can be obtained, according to considerations 
presented earlier, simply by letting ti — > oo in the general 
solution. Since Y is less than unity, it is evident that, for 
very large values of n, 


r C dY=l-Y, and 7„-^l. 

Using equations (19) and (25), and noting that Z=0 for the 
present case, it follows that 

^ f F on surface facing upstream 

(2 — F on surface facing downstream 

B=2{2-H) (29) 


1= 1/(2 -27) = 2/P 


(30) 


Equation (28) shows that the optimum surface has a dis- 
continuity in slope at X=\, and that both segments make 
a common angle with the chord plane. Equation (29) pro- 
vides the required relation between the base pressure param- 
eter B and the optimum trailing-edge thickness. Equation 
(30) determines the position of maximum thickness. As H 
varies from 0 to 1 equation (29) covers only the range of B 
from 2 to 4. Within this range the above equations apply, 
and the optimum profile is of the type illustrated in figure 
3. For the range of B from 0 to 2 the second possible solu- 
tion discussed earlier, namely, a wedge profile, represents 
the optimum section. For values of B greater than 4 the 
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double-wedge airfoil with a sharp trading edge is tbe optimum 
for a given thickness ratio. 

It is remarked that the above solution also can be ob- 
tained quite easily by solving equations (14) and (15) 
directly, instead of employing the limiting process. This 
direct method can be used to determine the optimum profile 
in those cases where the general solution in terms of n and a 
cannot readily be foimd. Such is the case when shock- 
expansion and similar higher-order theories are employed. 

CALCULATION OF PRESSURE DRAG COEFFICIENT 

Given structural criteria. — Since the dimensionless thick- 
ness distribution of an optimum profile is completely deter- 
mined by the base pressure parameter, it is to be expected 
that the quantity /3ci/(i5/c)^ also wdl depend only on B. 
From equation (1) and the definition of B it is seen that 


&Ci 2dc 

{ticy- 



h 

t 


Substituting P=2y'l^ and changing to the dimensionless 
variables X, Y, and H yields 


This equation gives the ratio or the pressure drag of an 
optimum profile to that of a sharp-trailing-edge, circular-arc 
profile having an equal value for I. If the pressure drag 
coefficient of a double-wedge profile (cg)a«, is used as a basis 
of comparison instead of a circular-arc profile, there results 
in a similar manner 

r, ky{l-[I/(t/2Y-n}+Bff(s/cy 
4(s/cy{ (u-f l)lI/(i/ 2 Y-^] 

It may be noted that the right side of equations (33) and 
(34) depend only on the base pressure parameter if the values 
of n and a are given. 

Given thickness ratio. — As noted earlier, as/i-^oo. 

From equation (27) it follows that 7/(</2)"“"-^0. By con- 
sidering equations (29), (30), and (31) there results 

Since [|SCtf/(//c)^]ca=16/3 and [|3Ctf/(f/c)^]tfu, = 4, it follows that 






The integral can be expressed in terms of 7/(f/2)"“‘" by 
noting that {dYIdKY^kyil — Y”-) and that 7/(t/2)”“‘'= 



F"7X 


There results for the pressure drag co- 


efficient of the optimum profile in linearized flow: 


/■ 2 1 

{t/cy " 


{sicy 


+BH 


(31) 


Inasmuch as 77, s/c, and II{tj2Y^'’ depend only on the base 
pressure parameter, the quantity on the left side of equation 
(31) also depends only on B for given values of n and c. 

It is of interest to compare the pressure drag of the theo- 
retically optimum profile with that of more conventional 
sharp-trailing-edge profiles. According to linear theory the 
drag coefficient of a biconvex circular-arc airfoil (q)^^ of 
thickness 7a is given by 


^{Ca)ca 16 

{tcjcy 3 


(32) 


A calculation of the value of the auxiliary integral for a 
circular-arc profile (lea) is readily made by substituting 
y=2(tcalc)x [1 — (a:/c)] in equation (2). It is found that 


7 ca I 


2^"in\y 


(2n-M)! 


(7a/2)"- 


and 

r Ti Ji- 

77^ = 9 -^=1-^ for2<5<4 (37) 

{Cd)m 2 16 4 — ^ ^ 

These are the same two equations that would be obtained by 
passing equations (33) and (34) to the limit as n-^ . When 

B<2, the optimum airfoil for a given thickness ratio, as 
previously discussed, is a wedge, for which PcJ(tlcy=l + B, 
Ctf/(cJca=3(l+5)/16, and cJ(c,),„=(l+5)/4. When B>4, 
the optimum is a double wedge, for which pcal(tlcy~4, 
cJ{Cd)ca=Sl4:, and cJ(Ctf)i„=]. 

If it is desired to compare the optimum profile with a 
corresponding sharp-trailing-edge profile on the basis of rela- 
tive 7 for a given Cg, rather than on the basis of relative 
for a given 7, then the foregoing calculations can be applied 
by making onij" minor modifications. As noted earlier, the 
thickness distribution of the optimum profile having maxi- 
mum 7 for a given is the same as that of the optimum 
profile having a minimum Ca for a given 7. By using the 
subscript s to denote a sharp-trailing-edge airfoil (e. g., bi- 
convex, or double-wedge), and no subscript to denote the 
optimum profile, the relation 

(L) 

\7s/c,^ = conBt. \^a/Z = const. 


By requiring that Ica=I where 7 is the value of the auxiliary 
integral for the optimum profile of thickness t and position 
of maximum thiclmess at s/c, then equations (31) and (32) 
can be divided to yield 


c, _ 7„^{l-[7/(f/2)"-^]}+,B77(s/c)^ 

(ca)ca ^ [7/«/2)«-1 


(33) 


can be deduced if it is remembered that 7 varies as the 
(n — <r) power of the thickness, and that the pressure drag in 
linearized theory varies as the square of the thickness. The 
above equation shows that in employing an optimum section 
the relative structural improvement that can be obtained 
for a given drag is related in a simple way to the relative 
drag reduction that can be obtained for a given structural 
requirement. 
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RESULTS FOR LINEARIZED FLOW AND DISCUSSION 

Significance and physical meaning of the base pressure 
parameter. — The determination of an optimum profile in 
linearized flow is greatly simplified by the fact that the 
dimensionless thickness distribution Y(x) depends only on 
the base pressure parameter B=—Pifil(tlc), and not on 
the individual values of P^, or tjc. 

Thus, although the Mach number, base pressui-e, and 
airfoil thickness ratio each indirectly affect the optimum 
airfoil profile, it is only necessary to Imow the value of B 
in order to determine the dimensionless thickness distri- 
bution {2y!t) of the optimum airfoil section. Knowledge of 
both I and B, of course, is sufficient to determine tjc as well 
as the dimensionless distribution of thickness. 

A simple physical interpretation of the base pressure 
parameter can be given if it is recalled that the basic means 
by which a thickened trailing edge reduces the over-all 
pressure drag is through a decrease in pressure foredrag at 
the expense of a smaller increase in hase drag. Thus the 
optimum dimensionless distribution of thickness must 
depend essentiallj^ on the ratio of base drag to pressure 
foredrag. The base drag for a given H is proportional to 
( — Pb)(tlc); whereas the pressure foredrag for a given V 
distribution is, according to linearized theory, proportional 
to (tlcyi-yjMJ-l. 

Hence, 

hase drag — Pt(t/c) _ —Pt,^|M„^—l _^ 

pressure foredrag (flcyi-yJM^^— 1 t/c 


or, in words, the base pressure parameter is proportional to 
the ratio of base drag to pressure foredrag. 

Condition under which optimum profile has a blunt trail- 
ing edge. — Prom equation (25) it is easy to deduce the con- 
dition under which the optimum airfoil will have a blunt 
trailing edge. The critical condition is obtained by setting 
H= 0. This determines a particular value of B, say B. 


■p 4nknin+2 — a-) 
{n—<r) (n+2) 


8 ... for 71= 1, <7 = 0 

6.283 . . . for 71=2, <r = 0 
5.609 . . . for 71=3, <7 = 0 
9.425 ... for n=2, <7 = 1 
6.730 . . . for 71=3, <7 = 1 
4 ... for 71= 00 , <7 finite 


(38) 


checks simple physical considerations. If the Mach number 
is infinite (for which B=0 even if a vacuum exists at the 
base), then this result checks the qualitative consideration 
of Saenger referred to in the introduction. 

A second limiting case that easily can be checked may be 
obtained by considering only the auxiliary conditions of 
given stiffness of sharp-trailing-edge profiles. The appro- 
priate results are obtained by setting <r=0 and H=0. From 
equation (26) it follows that Z=0. From equations (19) and 
(21) it is seen that the optimum sharp-trailing-edge profile 
is a doubly symmetric profile, each side of which is comprised 
of the arc of a parabola for ti=1, the arc of a trigonometric 
sine function for n=2, and the arc of an elliptic sine func- 
tion for 71=3; these are the results obtained previously by 
Drougge. 

Summary curves of the principal results. — In figure 4 the 
optimum dimensionless trailing-edge thickness H is plotted 
as a function of the base pressure parameter B. Each curve 
in this figure is obtained by substituting the indicated values 
of n and <7 in equation (25). It is to be remembered that 
the curve consisting of three straight-line segments, cor- 
responding to 71= CO and <7=finite, represents the auxiliary 
condition of a given thickness ratio. The other values of n 
and <7 represent the various structural criteria listed in the 
table presented earlier. 



Figure 4 —Optimum trailing-edge thickness for linearized supersonic flow. 


A lower value of B would correspond, for example, to a lower 
base drag, hence the physical significance of B can be stated 
quite simply; the optimum airfoil has a blunt trailing edge 
for B<CB; whereas it has a sharp trailing edge for B^B. 

Comparison with results of other investigations. — As. a 
partial check on the equations developed, several limiting 
cases can be obtained by specializing to particular values 
of n, a, and H. First, if the base pressure coefficient is zero, 
corresponding either to zero base drag or else infinite Mach 
number, then 5=0. From equation (25) it follows that 
H=l. In other words, the optimum profile for P},=Q has 
its maximum thickness at the trailing edge. If the Mach 
number is finite and the base drag zero, then this result 


The location of the optimum position of maximum thick- 
ness s/c, as determined by equations (24) and (26), is plotted 
in figure 5 as a function of B. The values of n and <7 used 
here are the same as in figure 4 . Comparing these two figures 
it can be seen that, as would be expected, the optimum posi- 
tion of maximum thickness moves steadily rearward as the 
optimum trailing-edge thickness is increased. 

Curves relating the value of I to the base pressure param- 
eter are shown in figure 6. These curves represent equa- 
tion (27). Since I is related to the optimum length of 
straight midsection through equation (13), the ordinate in 
this figure represents either of the two equal quantities, 
//(f/2)"“" or 7iZ/<7C. Figm'e 6, therefore, can also be used to 
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Figure 5.— Optimum position of maximum thickness for linearized supersonic flow. 



Figure 6 —Curves relating / and I to the base pressure parameter for linearized supersonic 

flow. 


determine lie in those cases where a is not zero. If <r is 
zero, then I is zero, as noted before. 

In figure 7 the two quantities cJ{Ca)ca and Cal(Ca)aw are 
plotted as a function of B for various values of n and <r. 
Depending on the value of B, it is apparent that the pressure 
drag of the optimum profile may be an 3 nvhere from a few 
percent to as much as 75 percent less than the pressure drag 
of an equivalent circular-arc sharp-trailing-edge airfoil. 
The structural criterion for which the greatest drag difference 
exists is that of a given bending strength of a thin-skin 
structure (n=2, o-=l). The curves of figure 7 (b) clearly 
illustrate the high drag of a double-wedge profile when it is 
compared to the optimum profile on the basis of a given 
structural requirement. These curves also illustrate that 
the relative drag reduction of the optimum airfoil for the 
condition of a given thickness ratio is much less than the 
corresponding reductions for the various conditions of given 
structural requirements. 

Method of determining an optimum profile from experi- 
mental base pressure data, — The experiments of reference 
8 have shown that the base pressure of airfoils in supersonic 




(a) Drag of optimum relative to biconvex circular-arc airfoils. 

(b) Drag of optimum relative to double-wedge airfoils. 

Figure 7. — Drag of optimum airfoils as compared to the drag of two different sharp-trailing- 
edge airfoils in linearized supersonic flow. Comparison made on the basis of equal values 
for the auxiliary integral. 

flow depends primarily on the trailing-edge thickness, 
Reynolds number, and type of boundary-layer flow. The 
base pressure generally does not depend significantlj^ on the 
shape of the airfoil profile upstream of the trailing edge. 
Figure 8 presents summary correlation curves (taken from 
reference 8) showing the dependence of base pressure on the 
parameters c/[A(i?e)'^®] and c/[h(ReY^^], which are approxi- 
mately proportional to the ratio of boundary-laj'er thickness 
to trailing-edge thickness for turbulent and laminar flow, 
respectivelj^ Since the optimum profile depends on the base 
pressure, which, in turn, depends on the trailing-edge thick- 
ness of the optimun profile, the value of B is not known 
initially. For this reason the process of determining an 
optimum profile from the experimental data of figure 8 
involves several steps: 

(1) For an arbitrarily selected value of t/c, B is computed 
as a function of H using the proper experimental value of 
base pressure for each H. 
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(a^ Turbulent, 



Figure 8. — Average values of base pressure from experiments of reference 8. 

(2) A plot of H versus B is superposed on figure 4. The 
point of intersection with the existing curve for the particular 
combination of n and a in question determines optimum 
values of H and B for the particular </c selected. 

(3) Knowing tjc, H, and B from (1) and (2), / is calculated 
from equation (27). 

(4) The above process is repeated for several values of tjc. 
Interpolation for the desired value of / then yields the values 
of tjc, H, and B of the optimum profile. The optimum 
value of B yields the optimum values of s/c and Z/c. Equa- 
tions (19) and (21) yield the basic shape of the curved 
portions of the desired airfoil. 

The results of applying steps (1) and (2) for a thickness 
ratio of 0.06, a Mach number of 3, and a turbulent boundary 
layer at i?e=10’', are shown in figiue 9. It is seen that the 
optimum trailing-edge thickness varies between 0.12f and 
0.67f for the different combinations of n and cr. The cor- 
responding pressure drag reduction compared to a biconvex 
airfoil having the same value for the auxihary integral varies 
between 6 and 29 percent, whereas compared to a double-wedge 


airfoil the corresponding pressure drag reduction varies 
between 1 and 63 percent. 

The effect of Mach number on the optimmn profile for 
i/c=0.06, Ti=l, (T=0, and turbulent flow at J?e=10^, is 
shown in figiue 10. For M„ = 5 an estimated value of 
p^jp„— 0.15 was. employed since experimental base pres- 
sure data are not yet available at this Mach niunber. For 
Mo>= “ and M„ = l, it is not necessary to know the base 
pressm’e to determine the optimum profile with linear theory. 
A large effect of Mach number on the optimum profile, 
particularly at Mach numbers above about 3, is evident 
from figiu-e 10. The effect of airfoil-thickness ratio on the 
geometry of the optimum profile also is large, as illustrated 
in figure 11. (For the case f/c=0.02 in this latter figure, it 
was necessary to extrapolate the experimental base-pressure 
curves of fig. 8 (a) in order to estimate the base pressure.) 
The trends illustrated in figm-es 10 and 11 can be explained 
from elementary physical consideration if it is recalled that 
B corresponds to the ratio of base drag to pressme foredrag. 
Thus, H approaches unity as approaches unity because 
the pressure foredrag in linear theory approaches infinity 
while the base drag remains finite. Moreover, H also 
approaches unity as Ma approaches infinity because the 
base drag, which is approximately proportional to 1/M^, 
becomes small compared to the pressm’e foredrag, which in 
linear theory becomes proportional to 1/M. By the same 
token, H approaches unity for very thick airfoils because the 
base drag, proprotional to tjc, again becomes small compared 



Figure 9. — Examples illustrating the effect of auxiliary. condition on the optimum profile; 
lin^rked flow, vertical scale expanded, Mco=S, f/c.=0.06, turbulent fl.ow at 
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Figure lO. — Examples illustrating the effect of Mach number on the optimum profile 
linearized flow, vertical scale expanded, ti = 1, <r=0, </c=0.06, turbulent flow at Re=10‘'. 


tjc^.oe 




to the pressure foredrag, which in linear theory is proportional 
to (</c)^. 

Reynolds number has an important effect on the optimum 
airfoil profile if the boundary layer is laminar. (See fig. 12.) 
This is because the base pressure depends markedly on 
Reynolds number for laminar flow. For turbulent flow the 
corresponding dependence is seen to be considerably less, 
and the optimum trailing-edge thickness is seen to be much 
less than for laminar flow. 


Re=!0 , Laminar 


H\ 

.75 


£ 





Re=!0‘, 

Laminar 







Re^!0% 

Turbulent 







Re=^/0\ 

Turbulent 
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Figure 12.— Examples illustrating the effect of Reynolds number and type of boundary-layer 
flow on the optimum profile; linearized flow, vertical scale expanded, n = l, ff=0, A/co =2, 
^/c*0.04. 


CONCLUDING REMARKS 

The general method presented for computing the profile 
shape having minimum pressure drag at zero lift has been 
developed for the auxiliary condition that 

^ = 7 f [2/70/2)1 dr 

c Jo 


is constant. For a given airfoil theory, the determination 
of an optimum profile under this condition involves the 
simultaneous solution of equations (11), (12), and (13), 
which are general in that the surface pressure coefficient 
F{y') and the parameters n and v are arbitrary. Such 
generality is useful since it allows either linear theory, 
second-order theory, or shock-expansion theory to be used 
in determining the optimum profile for a number of practical 
auxiliary conditions such as prescribed bending strength or 
given torsional stiffness. As an illustration of the method , a 
solution has been developed“in detail using linearized flow, 
that is, using the expression P=2y' j-y] M J^—1. In this 


Figure 1L— Examples illustrating the effect of airfoil thickness ratio on the optimum profile; 
linearized flow, vertical scale expanded, n=\, a=Q, Afoo=3, turbulent flow at Re=W. 
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simple case a complete solution in closed form is obtained 
for the thickness distribution of the optimum proffle. 

The principal result of the analysis for linear supersonic 
flow is that the dimensionless thickness distribution of the 
optimum profile depends only on the single parameter 

This parameter has been termed 
the base pressure parameter, and has a simple physical 
significance in that it is proportional to the ratio of base 
drag to pressure foredrag. The dependence of an optimum 
profile in linear flow on one parameter only enables summary 
cmves to be plotted showing all principal results as a func- 
tion of B (figs. 4, 5, 6, and 7). The optimum dimensionless 
trailing-edge thickness increases if either the base pressure 
is increased, the airfoil-thickness ratio is increased, or the 
Mach number is increased to very high values. 

At low supersonic Mach numbers the theoretical results 
obtained are questionable since the assumptions of linear- 
ized airfoil theory break down as the Mach number ap- 
proaches unity. The results can be applied safely only 
to cases where linear theory satisfactorily predicts the 
pressure foredrag. Although at high supersonic Mach num- 
bers the results obtained under the assumption of linearized 
flow also would not be expected a priori to be of quantita- 
tive value, they predict, nevertheless, the correct result 
that the optimum trailing-edge thickness for infinite Mach 
number is equal to the maximum airfoil thickness. In 
view of this exact agreement in the extreme case, it is con- 
jectured that the linear theory fortuitously may provide 
a reasonable estimate of the optimum trailing-edge thiclmess 
for any supersonic Mach number not close to unity. As 
regards the optimum profile shape forward of the base, 
however, such fortuitous conditions cannot be expected. 


since the hnearized approximation at high Mach numbers 
overestimates the suction forces and underestimates the 
positive pressure forces. This causes the calculated opti- 
mum profile to have too large a leading-edge angle, a posi- 
tion of maximum thickness too far forward, and too small 
an inclination of the surface behind the position of maxi- 
mum thickness. (In reference 9 some calculations using, 
second-order theory are presented which illustrate this 
effect on the optimum sharp-trafling-edge profile for the 
auxiliary condition of a given thickness ratio.) 

Because the optimum profile, by definition, has the 
least pressure drag possible under given conditions, small 
changes in profile shape would result in second-order changes 
in drag. This allows some flexibility in modifying the 
theoretically optimum profile to more closely suit individual 
design requirements, and means that it is not important 
to rigorously adhere to the exact parabolic, trigonometric- 
sine, or elliptic-sine contour (provided, of course, that the 
end points of the modified contour are located approximately 
in the optimum positions). It is important to adhere 
reasonably close to the calculated optimum trailing-edge 
thickness, since this quantity can greatly affect the drag. 
In particular, a trailing-edge thickness considerably greater 
than the optimum should not be used. Excessive trailing- 
edge thickness at low and moderate supersonic Mach 
numbers can result in an excessive increase in drag. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., October 3, 1960. 
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APPENDIX 


NOTE REGARDING DISCONTINUANCE CHANGE IN SLOPE 
OP AN OPTIMUM PROFILE 

In the general analysis of optimum profiles for a given 

dP 

thickness ratio it was foimd that the function y'^ necessar- 

ily was constant along each straight segment of the profile 
(equation (14)). No information was obtained, however, 
about the relative value of this constant for the two seg- 
ments. The required information can readily be obtained 
by considering the change in drag due to a change only in 
position of maximum thickness, that is, a change in slope of 
both straight surfaces illustrated in figure 3 with no change in 
t or h. Using subscripts u and d to denote surfaces facing 
upstream and downstream, respectively, it follows that 


— hic^u “b 5(l'd)(/ 


Since y\—tls and y'i=(t—h)l(c—s), the minimizing of Ca 
requires that 


(dy'A (c-s)^ 




2 

y d 


(i^)j 


bP 

Thus 2 /'^ must 
" 5 ?/ 


be continuous at the corner — a result 


which was used without proof in the general analysis. It 
may be noted that an alternate proof of this result can be 
obtained in an extremely easy way from the following 
known result of the calculus of variations: The Weierstrass 
E-function is continuous at the point of discontinuity on a 


boundary. The E-function in the present case is y'^ 

For auxiliary conditions other than a given thickness ratio 
it was tacitly assumed in the analysis that the optimum 
surface everywhere had a continuous slope. This assumption 


also requu-es some justification. It is shown in the calculus 
of variations that at all points of free variation it is necessary 

for to be continuous. For the present problem this means 
bP 

that P+y' must be continuous at all such points. Ac- 

cording to linear theory, P+y' hence, within the 

scope of linear theory, the surface slope y' is continuous at all 
points of free variation. For shock-expansion theory 
bP . 

is positive, and a corner would cause a discontinuous de- 
bP 

crease in P, y', and P+y' hence, also within thescopeof 

shock-expansion theory the surface slope of the optimum pro- 
file is continuous at all points of free variation. This justifies 
the assumption of continuous slope employed in the general 
analysis for auxiliary conditions other than a given thickness 
ratio. 
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